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Abstract 

Results are presented from an application of H^ control design methodology to a 
centralized integrated flight /propulsion control (IFPC) system design for a supersonic 
Short Take-Off and Vertical Landing (STOVL) fighter aircraft in transition flight. The 
emphasis is on formulating the H^ control design problem such that the resulting controller 
provides robustness to modelling uncertainties and model parameter variations with flight 
condition. Experience gained from a preliminary H^ based IFPC design study performed 
earlier is used as the basis to formulate the robust H^ control design problem and improve 
upon the previous design. Detailed evaluation results are presented for a reduced order 
controller obtained from the improved H^ control design showing that the control design 
meets the specified nominal performance objectives as well as provides stability robustness 
for variations in plant system dynamics with changes in aircraft trim speed within the 
transition flight envelope. A controller scheduling technique which accounts for changes in 
plant control effectiveness with variation in trim conditions is developed and off— design 
model performance results are presented with the scheduled controller. 
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Introduction 

The trend in future military fighter/tactical aircraft design is towards aircraft with 
new/enhanced maneuver capabilities such as Short Take-Off and Vertical Landing 
(STOVL) and high angle of attack performance. An integrated flight /propulsion control 
(IFPC) system is required in order to obtain these enhanced capabilities with reasonable 
pilot workload. An integrated approach to control design is then necessary to achieve an 
effective IFPC system. Such a design approach is currently being developed at NASA 
Lewis Research Center under an in-house research effort. This methodology is referred to 
as IMP AC - Integrated Methodology for Propulsion and Airframe Control [1]. The 
significant features of the IMP AC methodology are that it consists of first designing a 
centralized controller considering the airframe and propulsion systems as one integrated 
system and then partitioning the centralized controller into decentralized subsystem 
controllers for state-of-the-art IFPC implementation. 

The major issue related to the centralized controller design portion of IMPAC is the 
choice of the control synthesis technique that "best" suits the IFPC objectives. Not only 
should the synthesis technique provide for means of formulating the design criteria for the 
centralized control design such that it adequately reflects the performance specifications of 
the "total" system, i.e., the airframe integrated with the propulsion system, but it should 
also result in controllers of reasonable complexity with guaranteed performance and 
robustness characteristics. Recent advances in control theory [2,3] and computational 
algorithms to solve for optimal control laws [4] have made this theory a viable 
candidate to be applied to complex multivariable control design problems. A preliminary 
investigation of the applicability of H control theory to the centralized feedback controller 

CD 

design portion of the IMPAC approach was conducted earlier via an example IFPC design 
study [5]. The results reported in Ref. [5] are quite encouraging in that they demonstrate 
that control theory has the promise to provide the framework to meet the requirements 
of a centralized IFPC design. However, the control design of Ref. [5] was a preliminary 
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design and detailed evaluation of that design identified various areas in which the 
controller performance needed to be improved. The objectives of this paper are to develop 
insight into formulating a robust control problem within the control design framework 
and to improve upon the control design of Ref. [5]. Towards this goal results are presented 

from a H based IFPC redesign for the linear model of a STOVL aircraft in transition 

00 

flight considered earlier in Ref. [5]. Also presented in this paper are the closed-loop 
performance evaluation results for off— design plant models using controller scheduling 
based on a scheme that exploits the robustness properties of the controller that are 
built into the nominal control design by appropriate problem formulation. 

The paper is organized as follows. The vehicle models to be used for control design 
and evaluation are first discussed. The control design is then presented along with some 
discussion of the formulation of the IFPC design objectives within the framework of the 
control problem. The emphasis is on formulating the problem such that the resulting 
controller is robust to modelling uncertainties and parameter variations with changes in 
flight condition. Evaluation results are presented for a reduced order approximation of the 
improved H controller and improvements over the controller design of Ref. [5] are 

CD 

demonstrated. A "simple" controller scheduling scheme that accounts for changes in plant 
control effectiveness for varying trim conditions is then discussed and evaluation results are 
presented for off— design models. 

Vehicle Model 

The vehicle considered in this study is representative of the delta winged E— 7D 
supersonic STOVL airframe powered by a high bypass turbofan engine [6]. The aircraft is 
equipped with the following controls: ejectors to provide propulsive lift at low speeds and 
hover; a 2D— CD vectoring aft nozzle with afterburner for supersonic flight; a vectoring 
ventral nozzle for pitch control and lift augmentation during transition; and jet reaction 
control systems (RCS) for pitch, roll and yaw control during transition and hover. A 
schematic diagram of the aircraft with relative location of the various control effectors 


4 


mentioned above is shown in Fig. 1. Engine compressor bleed flow is used for the RCS 
thrusters and the mixed engine flow is used as the primary ejector flow. Detailed ducting 
diagrams of the engine and discussion of the ejector STOVL concept are available in Ref. 
[ 6 ]- 

Currently, two separate computer simulations, one for the airframe and one for the 
propulsion system, are being used to assess performance capabilities of the aircraft and to 
generate open— loop linear models for control design [6]. The procedure for generating 
integrated airframe and propulsion models for control design and evaluation from the two 
separate simulations is discussed in Ref. [5]. The integrated linear design model used in this 
study is of the form 


x = Ax + Bu ; y = Cx + Du (1) 

where the state vector is 

x = [u,v,w,p,q,r,^,0,N2,N25,T41,T3,P6] T (2) 

with 

u = Axial Velocity, ft/s 

v = Lateral Velocity, ft/s 

w = Vertical Velocity, ft /s 

p = Roll Rate, rad/s 

q = Pitch Rate, rad/s 

r = Yaw Rate, rad/s 

<t> = Roll Attitude, rad 

6 = Pitch Attitude, rad 

N2 = Engine Fan Speed, rpm 

N25 = High Pressure Compressor Speed, rpm 

T41 = High Pressure Turbine Inlet Temp., °R 

T3 = Combustor Inlet Temp., °R 

P6 = Tailpipe Entrance Total Pressure, psi, 

the control inputs are 

u= [MMr,AQR,AYR,ARR,WF,A8,ETA,A78,ANG8,ANG79] T (3) 

with 


& = Elevator Deflection, deg 

& = Aileron Deflection, deg 

Si = Rudder Deflection, deg 
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AQR = Pitch RCS Area, in 2 

AYR = Yaw RCS Area, in 2 

ARR = Roll RCS Area, in 2 
WF = Fuel Flow Rate, lbm/hr 

2 

A8 = Aft Nozzle Throat Area, in 
ETA = Ejector Butterfly Angle, deg 

2 

A78 = Ventral Nozzle Area, in 

ANG8 = Aft Nozzle Vectoring Angle, deg 

ANG79= Ventral Nozzle Vectoring Angle, deg, 

and the outputs are 


y = [V,V,«, q> 7,#,P,flr,N2,fl T 


( 4 ) 


with 

V = True Airspeed, ft/s 

2 

V = Total Acceleration, ft/s 

7 = Longitudinal Flight Path Angle, deg 

fa = Sideslip Angle, deg 

/? = Rate of change of Sideslip Angle, deg/s. 

The other outputs are as discussed under state description except that the angular 
positions and rates are in degrees. 

Some discussion of the choice of control inputs u for the linear design model is 
relevant here. The choice of outputs y is discussed later in the paper. The E— 7D aircraft is 
equipped with left and right elevons on the trailing edge of the delta wing. Collective 
deflection of the elevons provides the classical elevator pitch control while differential use 
of the elevons provides the aileron roll control. So the elevator (£e) and aileron (6a) along 
with the rudder (tfr) are used as the airframe control inputs in the design model. Only 3 
RCS areas, AQR, AYR and ARR, are used as RCS control inputs in the linear design 
model whereas the full nonlinear model has 5 controlled RCS areas. The reasons for this are 
that the nose pitch RCS and the two yaw RCS thrusters provide thrust in only one 
direction as shown in Fig. 1, and the wing tip roll RCS thrusters are to be used 
differentially for roll control and collectively for pitch control. For instance yaw RCS 
thrusters provide only forward thrust, so left yaw RCS is used for right yaw and right yaw 
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RCS is used for left yaw in the nonlinear model. Using both left and right yaw RCS areas 
as control inputs in the design model can result in a control design that uses the two areas 
differentially to enhance yaw control which is inconsistent with the actual implementation. 
Similar reasons apply for pitch and roll RCS area selections. An RCS distribution logic 
that will distribute the three design model RCS commanded areas to the five actual areas 
in the nonlinear control implementation is shown in Fig. 2. Since the nose pitch RCS 
thruster only provides positive (pitch up) pitching moment, a negative AQR command is 
distributed among the left and right wing tip RCS thrusters to generate the commanded 
negative pitching moment taking into account the relative pitch control effectiveness of the 
nose and wing tip RCS thrusters. 

Fig. 2 is presented here to clarify the discussion on choice of RCS areas for linear 
design model although the logic of Fig. 2 is not used to evaluate the control designs 
presented in this paper. Note that there is an absolute nonlinearity in the relationship from 
commanded RCS areas to compressor bleed flow demand in that although the RCS area 
may be positive or negative depending on the desired direction of RCS thrust, the 
compressor bleed flow demand (WB3) to generate the thrust is always positive. For a 
linear model trimmed about zero RCS areas, this relationship is of the form 
WB3i Kj-| AiR| where AVB3i is the demanded bleed flow due to AiR command with i 
representing Q, Y or R for pitch, yaw or roll RCS, respectively, K . is an appropriate 
constant and | • | represents absolute value. This nonlinearity was not taken into account 
in the linear control design of Ref. [5] and the resulting controller was shown to lead to 
unacceptable performance when evaluated with this nonlinearity. In this paper, the critical 
absolute nonlinearity from RCS area commands to WB3 demand is accounted for in the 
linear control design process by a proper formulation of the control design problem within 
the framework. The details of this formulation are discussed later in the paper. 

The control inputs, WF, A8, ETA, A78, ANG8 and ANG79 in the vector u, are the 
propulsion system inputs. The ejector butterfly valve angle (ETA) controls the engine 
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airflow to the ejectors, thus providing a means of controlling ejector thrust. There are 
separate control valves for the left and right ejectors, however the two valve angles are set 
to be equal in the aircraft simulation because no test data is available on the differential 
use of the ejectors for roll control. Therefore only one butterfly valve angle is used as the 
control input in the design model. The other five propulsion system control inputs in the 
design model were just as defined in the full cycle-deck engine simulation [6]. 

The flight phase considered in this study is the decelerating transition during 
approach to hover landing. This phase presents a challenging control design problem 
because the control of the aircraft is transitioning from aerodynamically generated forces 
and moments to those generated by the propulsion system. For this study, three linear 
integrated models were obtained corresponding to steady-state level flight at trim speeds 
of V Q = 60, 80 and 100 Knots with trim flight path angle 7 Q = —3 deg for all three models. 
The trim strategy used to generate the linear models is as follows : i) The aircraft pitch 
attitude was set at # Q = 7 deg to provide adequate visibility during landing, and the elevon 
settings were chosen to correspond to elevator deflection & Q = 20 deg so that adequate 
elevator control authority is available for maneuvers during transition; ii) The aft nozzle 
vectoring angle was set to ANG8 Q = 0 so that the total vectoring authority is available for 
pitch control, and the ventral nozzle vectoring angle was set to ANG79 Q = 64 deg to ensure 
that the aircraft can be trimmed at the low speed (60 Knots) with small but nonzero aft 
nozzle thrust, hence nonzero throat area. This trim strategy leaves adequate A8 and 
ANG79 control authority available for these controls to be used for active speed and pitch 
control in transition; iii) The aircraft is then trimmed by using the three thrusts — aft 
nozzle, ventral nozzle and the combined ejectors, as the trim controls. The resulting trim 
values of A8, ETA, and A78 for the three trim speeds are listed in Table 1. Also listed in 
Table 1 are the trim values of fuel flow (WF) and engine fan speed (N2) which define the 
engine operating point. Note that the quantities shown in Table 1 are a percentage of the 
corresponding trim values for the 80 Knot model. As seen from Table 1, thrust is 
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transferred from aft nozzle to ventral nozzle as the aircraft slows down and the engine 
operating point is moved up to a higher gross thrust level (increased WF and N2) to 
generate the propulsive lift needed to compensate for the loss in aerodynamic lift. 

Open loop frequency— domain and time— domain analyses of the three linear models 
indicated that the 80 Knot model provides a "good average" of the three models in terms of 
input/output response behavior. Thus the 80 Knot model was used for control design with 
the 60 Knot and 100 Knot models being used to evaluate the stability and performance 
robustness of the nominal design, develop controller scheduling and evaluate off— design 
performance. The eigenvalues of the linear models for the three flight conditions are listed 
in Table 2 and the airframe modes are identified in terms of their "classical" interpretation 
[7]. As seen from Table 2, the aircraft is unstable in pitch with the short period (SP) mode 
becoming less unstable as the trim speed decreases. Also the dutch roll mode damping is 
very low and it decreases with trim speed with the mode going unstable at 60 Knots. The 
engine temperature and pressure modes are very fast compared to the engine rotor 
dynamics and aircraft modes. Ideally these fast modes will not be considered in the linear 
design model in order to reduce the complications in the control synthesis procedure. 
However, the approach taken in this study is to include these modes in the design model 
and reduce the complexity of the controller through order— reduction after initial controller 
synthesis. 

Control Design 

Design Methodology 

The IFPC design problem discussed earlier was formulated as a command tracking, 

disturbance rejection problem within the framework of the general mixed sensitivity H 

control problem [8]. The detailed block diagram for the H formulation of the IFPC 

00 

feedback control design is shown in Fig. 3. The three transfer functions that are of interest 
for such a problem are the sensitivity function S(s), the complementary sensitivity function 
T(s), and the control transmission function C(s). These represent the closed-loop transfers 
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from the reference commands and disturbances to tracking errors, controlled variables and 
commanded control inputs, respectively. In order to influence both the low-frequency and 


high-frequency properties of the closed-loop system it is desirable to find a controller K(s) 
which minimizes a weighted norm of a combination of these three transfer functions, i.e. 

rw s a«)-s(jw)i 


min||H(jw)|| (B with H(jw)= 


W T (jw)*T(jw) 

W c (j«)-C(j«) 


Note that ||H(jw)|| is the highest value over all frequencies u of the maximum singular 
value of H(jw), i.e. 

PCMI = max RH(jw)]) 

u 

The weighting functions Wg, Wrp, and W^jw) are the "knobs" used by the control 
designer to "tune" the controller K(s) such that the design objectives are met. For instance 
choosing Wg to be large at low frequency ensures good command tracking performance and 
choosing W,p to be large at high frequencies ensures robustness to high frequency 
unmodelled dynamics. are chosen to ensure that control actuation bandwidths and 
control rate and deflections are practically achievable. 

The tracking formulation of Fig. 3 allows for feedback of plant measurements 
other than just tracking errors as inputs to the controller. This formulation then allows the 
simultaneous design of inner loop plant augmentation (stability or response shaping) and 
command tracking system. Such plant augmentation is an integral part of flight control 
design since the overall objective is to design a system for desired piloted handling qualities 
and not just an automatic command tracking system. Also, it has been shown in Ref. [9] 
that if a tracking problem is formulated within the framework as purely a 
servo-mechanism problem, i.e. with controller inputs being just the tracking errors, then 
the controller will be such that its transmission zeros cancel the stable poles of the 
design plant thus resulting in a closed— loop system that will have almost no robustness to 
plant modeling uncertainties. Allowing for feedback augmentation within the H control 
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formulation overcomes this problem of lack of robustness as demonstrated in Ref. [5]. 

Design Specifications 

The vectors u and y in Fig. 3 are the integrated design model inputs and outputs, 
respectively, as discussed in the previous section. The controlled variables z were selected 
to be 

z = [Vv,Qv,7,Pv,/?,N2] T (5) 

with Vv=V+0.1V, Qv=q+O.30, Pv=p+0.1<^ and the others as discussed under plant 
outputs. This blending of controlled variables was chosen to provide the response types 
that are desirable for good handling qualities [10,11] in transition flight. The choice of Vv 
corresponds to designing an acceleration (deceleration) command system with velocity 
hold, and the choice of Qv and Pv correspond to designing a rate command— attitude hold 
system. The break frequencies for switching from rate to attitude command for the case of 
Qv and Pv, and from acceleration to velocity command for the case of Vv, were chosen 
based on open-loop control effectiveness studies. For instance, the elevator (<fe) is effective 
in pitch rate control in the frequency range of 0.3 rads/s to 10 rads/s and is effective in 
pitch attitude ( 6 ) control for frequencies below 0.3 rads/s. The choice of 7 in z provides for 
flight path angle control and the choice of N2 provides for tracking the fan speed 
commands generated by the engine operating schedule logic. 

Designing the feedback controller K(s) to provide decoupled command tracking of 
the individual elements of z will result in a system that provides independent control of 
acceleration, pitch, flight path angle, roll and sideslip from the various pilot control 
effectors such as stick, throttle and rudder pedals etc., thus reducing pilot workload, and 
also control of the propulsion system operating point (N2) independent of the aircraft 
motion. Independent control of roll (Pv) and sideslip angle (/?) will result in a control 
system that provides automatic turn coordination thus further reducing pilot workload. 

The choice of outputs y, in (4), is consistent with the various quantities used for 
feedback in classical aircraft control augmentation [7]. Rate feedback provides for improved 
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damping while position (angles) feedback improves natural frequency of the appropriate 
mode. 

Based on the experience gained from the control design study of Ref. [5], an 
important criterion for control design is that the closed— loop system be robust to the 
absolute nonlinearity associated with the bleed flow demand from the RCS area commands — 
see earlier discussion. Also, the control design of Ref. [5] was such that it did not further 
stabilize the dutch roll mode, i.e. the open— loop dutch roll mode was also a mode of the 
closed— loop system with the designed controller. This resulted in an unstable closed— loop 
system for the off— design integrated model at 60 Knots due to the dutch roll mode being 
unstable at that flight condition. For this redesign effort, stabilization of the dutch roll 
mode to guarantee stable closed loop system over the transition flight conditions being 
considered was included as a design specification. 

Control Design 

The design plant inputs and controlled outputs were normalized by maximum 
allowable deflections (u ) and maximum commanded values to be tracked (z ) 

v max' v C in ax' 

respectively. The u were chosen to be reasonable deviations from the nominal (trim) 

max 

values such that the total control deflection limits (as incorporated in the actuator models 
in the nonlinear model) will not be exceeded. Some of the u values corresponding to the 

max 

propulsion system controls were further reduced to ensure that the total control usage will 
be within the limits imposed by safety requirements. For example, the allowable engine 
fuel flow for a given operating point is limited by the engine acceleration/deceleration 
schedule. The z were chosen based on handling qualities control requirements and 
open— loop control effectiveness analysis of the design plant to ensure that each element of z 
can be commanded individually to its maximum value within its frequency range of 
interest without exceeding u value for any of the control inputs. The singular values of 

max 

the scaled design plant for the 6 controlled variables defined in eqn. (5) are shown in Fig. 4. 
The fact that the minimum singular value in Fig. 6 is less than 1 implies that there exist 
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combinations of numerical values of commands z such that although each element of z is 

c c 

less than its maximum value, the combined commands z cannot be tracked without 

c 

exceeding the control limits u for some control input. 

max 

The sensitivity weights Wg and the complementary sensitivity weights Wrp for each 
of the controlled variables were chosen to be first — order, to provide adequate frequency 
response shaping without overly increasing the resulting controller order. The Wg zero and 
pole for each controlled variable were chosen to result in a low frequency gain of 1000, gain 
crossover frequency of 3 to 4 times the control bandwidth desired for good handling 
qualities for the controlled variable of interest, and a high frequency gain of 0.1. Such a 
choice of Wg reflects the desire to synthesize a sensitivity function which give good 
steady— state tracking performance in the presence of disturbances and low frequency 
modelling errors, good tracking performance up to the desired bandwidth of control and 
reduced emphasis On tracking at high frequencies where there are significant modelling 
errors and uncertainties. The Wrj, were chosen to obtain a low frequency gain of 0, gain 
crossover frequency of about 1.2 times the corresponding Wg gain crossover frequency, and 
a high frequency gain of 100. This choice of W,j, ensures that the plant outputs are not 
penalized at low frequencies where command tracking is to be emphasized while at high 
frequencies the plant outputs are penalized heavily to provide controller gain attenuation 
for robustness to unmodelled dynamics. 

As part of the control weighting Wq, both the control inputs and control rates were 
weighted with W in Fig. 3 chosen to be the inverse of u and W • chosen to be inverse of 

max U 

maximum control rate for each control input. Since using full order actuator models for 
each control input would have resulted in a very high— order controller, first order actuator 
approximations were used in the control design. Describing function analysis [12] of the full 
order actuators was first performed to determine the degradation in actuator bandwidth 
due to rate limiting when control commands corresponding to u are used at all 

max 

frequencies. The worst-case rate-limited actuator bandwidth was then used as the 



13 


bandwidth for first-order design actuators. For example, the pitch RCS area actuator 

bandwidth is 20 rad/s but the rate limit is 3.0 in /s which results in a worst case 

2 

bandwidth of 6.0 rad/s for AQR of 0.7 in . Using such an approximation for the 

max 

actuators in the design guarantees performance robustness in the presence of actuator rate 
limiting. 

The bleed flow demand from the RCS acts as a disturbance on the engine operating 
point. In the propulsion system simulation and the linear model generated from this 
simulation, bleed flow appears as an external disturbance. However, when the airframe and 
propulsion models are integrated there is no explicit dependence on bleed flow since the 
bleed flow is generated from' the RCS area commands which are the inputs to the 
integrated model. As pointed out earlier, the linear design model does not account for the 
absolute nonlinearity from the RCS area commands to bleed flow demand, so the effect of 
the RCS areas on the engine states as it appears in the linear design model is erroneous. 
The linear design model was modified such that the elements of the B matrix from the RCS 
areas to engine states are zero, i.e. the RCS areas do not affect the engine operating point, 
and bleed flow was added as an external disturbance affecting the engine dynamics as in 
the linear engine model. For control synthesis, the bleed flow disturbance was modelled 
as a zero— mean Gaussian white noise process with an intensity of 7 lbm/s, which is the 
maximum possible RCS bleed flow, and was considered to be a measurement available for 
feedback. Although bleed flow is not measurable, it can be estimated as a function of RCS 
areas and ambient and bleed duct pressure conditions or the bleed flow feedback can be 
implemented as nonlinear feedback of the RCS area commands. Since the RCS is used for 
control of aircraft angular rates, the bleed flow disturbance was filtered by a first order 
filter with a bandwidth of 3.5 rads/s which is representative of desired angular rate control 
bandwidths. 

A procedure for providing robustness to plant model variations within the 
framework of the control problem is to use internal noise models as disturbances which 
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mimic the effects of model variations [13]. Such an approach was used in the current design 
study to provide enhanced lateral/directional stability by appropriate damping of the 
dutch roll mode and robustness to changes in the lateral/directional dynamics with 
decreasing trim airspeed. Side force, rolling moment and yawing moment changes were 
represented in the form of zero— mean white noise Gaussian disturbances in lateral 
acceleration (v), roll acceleration (p), and yaw acceleration (r) with intensities of 0.5 ft/s^, 
0.1 rad/s , and 0.1 rad/s respectively. The overall control synthesis problem was then 
formulated as that of command tracking with bleed flow and lateral/directional 
disturbance rejection. 

The design plant as discussed above is of 38th order consisting of the 13th order ' 
integrated airframe/propulsion system design model, first order actuators for the 12 control 
inputs, first order sensitivity and complementary sensitivity weights for the 6 controlled 
variables, and the first order filter for the bleed flow disturbance. So the H controller 

0D 

using the algorithm of Ref. [2] will be of 38th order. The control design results with this 
design plant are shown in Fig. 5 in terms of the closed loop maximum singular values of the 
combined weighted functions <F[H(ja;)], weighted errors <r[Wge(ju;)], weighted controlled 

outputs a[ Wrpz(jw)], weighted controls a[W u u(jw)], and weighted control rates a[W^u(jw)] 
with the commands z £ as inputs. ||H(jfa/)|| s 10 is achieved for this controller as seen from 
Fig. 5. In general a control design with ||H(jt«/)|| < 1 ensures that all the design 
specifications that are formulated through the various weightings will be met. However, 
this was not the case in the present design, because as stated earlier the minimum singular 
value of the scaled design plant itself was much less than 1. The fact that control efforts 
greater than u will be required to track some combinations of commands is evident from 

max 

the large (>1) maximum singular values of weighted controls at low frequencies. The fact 
that maximum singular value of weighted errors is greater than 1 at low frequencies 
indicates the performance/control trade-off made in the H minimization procedure. The 
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fact that the maximum singular values of the weighted controlled outputs and control rates 
are less than 1 over all frequencies indicates that the controller provides adequate "loop 
gain" roll-off for the closed— loop system to be robust to unmodelled high frequency 
dynamics and that the control rate limits will not be exceeded for any combinations of 
commanded variables. Further evaluation results are presented in the following section 
with reference to a reduced— order controller obtained from the H controller. 

QD 

Controller Reduction and Evaluation 

The controller order was reduced from 38 to 14 using a combination of modal 
residualization [14] and balanced realization [15] reduction techniques. The order was first 
reduced from 38 to 22 by modal residualization of the high frequency controller modes. The 
22 n( ^ order controller was reduced to 16^ order by a balanced realization approximation 
which was then reduced to 14^ order again by residualization of the high frequency modes. 
The eigenvalues of the final reduced order controller were bounded by | A | <10 thus 
indicating that digital implementation of the controller can be achieved with reasonable 
sampling rates. The performance with the full order and reduced order controllers is 
compared in Fig. 6 in terms of the maximum and minimum singular values of the 
closed— loop tracking system, a[T(jw)] and a[T(jw)] with z(s)=T(s) z c (s). Clearly, there is 
no significant loss in tracking performance with the reduced order controller. Note that 
other modern control reduction techniques such as frequency— weighted balanced realization 
[16] were not used because it was not clear how to apply these techniques for a controller 
structure which includes external disturbances as inputs (bleed flow in this particular case). 

Extensive frequency— domain and time— domain analyses were performed to evaluate 
the closed— loop performance and robustness with the linear design model. All these 
analyses indicated that the reduced-order controller provides decoupled command tracking 
with desired tracking bandwidths and reasonable control actuation requirements. Note that 
although the open— loop design model is decoupled in the longitudinal and 
lateral/directional dynamics, coupling between these axes is introduced through the use of 
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RCS control. For instance using the yaw RCS will result in change in velocity because the 
yaw RCS generates axial force. So the fact that the control design achieves decoupled 
response in longitudinal and lateral/directional axes while using RCS for active control is 
significant. 

The performance of the control design in the presence of the absolute nonlinearity 
from commanded RCS areas to compressor bleed flow demand was evaluated by including 
these nonlinear effects in the linear 80 Knot design model. An example result is shown in 
Fig. 7 in terms of the closed— loop pitch attitude and engine fan speed response to a pitch 
variable (Qv) command. The Qv command was chosen to correspond to a transient pitch 
rate command and a steady— state pitch attitude command. Also shown in Fig. 7 is the 
corresponding closed— loop system response with the control design of Ref. [5], The current 
design provides the desired improvement over the control design of Ref. [5] in that it 
maintains good pitch tracking and fan speed regulation performance in the presence of the 
absolute bleed flow nonlinearity. Thus the formulation of the control problem to reject 
bleed flow disturbance with the modified design plant, the modification being the removal 
of the RCS area effect on the engine dynamics, was successful in providing performance 
robustness to the bleed flow absolute nonlinearity. 

Another important design criterion discussed earlier was to provide stability 
robustness to changes in the lateral/directional dynamics specially with reference to the 
low damping of the dutch roll mode. For the current control design, the open-loop dutch 
roll mode was no longer a closed-loop mode with the 80 Knot design plant. Moreover, all 
the closed— loop modes were well damped with the smallest damping ratio being 0.4. 
Detailed stability robustness evaluation of the control design to variations in the plant 
system A matrix due to changes in trim speed were conducted using structured singular 
value robustness tests [17]. The procedure for this stability robustness evaluation is 
discussed in the following. 

First the variations in the airframe portion of the plant A matrix with changes in 
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trim speed were identified. With the 80 Knot model A matrix as the nominal, thirteen 
elements of the airframe portion of the A matrix showed a change of 25 % or greater in 
going from 80 Knots to 60 Knots or 80 Knots to 100 Knots. These elements correspond to 
change in axial (X-axis) force due to pitch rate, side force (Y-axis) due to roll rate and 
yaw rate, vertical force (Z-axis) due to axial speed (u), vertical speed (w) and pitch rate, 
rolling moment due to roll rate and yaw rate, pitching moment due to axial speed and 
pitch rate, and yawing moment due to side velocity (v), roll rate and yaw rate. Also the 

percentage change in the A matrix elements was approximately symmetrical about the 

nominal A matrix. For a given trim speed, Vrp, these thirteen elements of the A matrix 
were then modelled using the multiplicative uncertainty form as 

Ajj ~ Ajj (l+Ky j AV) (6) 

where A°j is the i,j element of the 80 Knot model A matrix, AV=V,j— 80, with units in 
Knots, and Ky^ were chosen to be the average of the change in the A(i,j) element from the 
80 Knot to 60 Knot and 100 Knot models. With such an uncertainty model, the 
perturbation A(s) for robustness analysis is of the form A(s) = AV-I where I is an identity 
matrix of dimension 13, and the structured singular value stability robustness measure p is 
given by [18] 

H = max [p(M(jw)] (7) 

u: 

where p[A] is the spectral radius of matrix A (maximum absolute eigenvalue of A) and 
M(s) is the interconnection matrix for robustness analysis which represents the nominal 
closed— loop system taking into account the structure of A(s). A detailed discussion of 
creating the interconnection matrix for this particular study is beyond the scope of this 
paper. The reader is referred to Ref. [17] for a theoretical discussion of the subject of 
robustness analysis for structured uncertainties and to Ref. [4] for a detailed exposition on 
creating the interconnection matrix for various forms of uncertainty models. The 
structured singular value for A matrix uncertainty as modelled by Eq.(5) is shown in Fig. 
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8. Also shown in Fig. 8 is the corresponding result with the control design of Ref. [5]. From 
Fig. 8 , n = .02 for the current control design which implies that for the A matrix 
uncertainties as modelled, the closed— loop system will remain stable for velocity 
perturbations of up to AV = i = 50 Knots. For the design of Ref. [5], \i = 0.065 which 
implies that the closed loop system with that control design would have been guaranteed 
only for velocity perturbations up to 15 Knots. In fact, for AV = 15 Knots as modelled, the 
open loop A matrix is such that the dutch roll mode is just unstable and the control design 
of Ref. [5] was found to be unstable for this perturbation. The closed-loop system with the 
current control design was stable for the full off— design 60 Knot and 100 Knot plant models 
thus indicating that the use of internal noise model in the H problem formulation was 
successful in building in robustness to variations in the lateral/directional dynamics. 

Since the design plant is an integrated flight propulsion system, robustness of the 
design closed— loop system to changes in the propulsion system dynamics is also an 
important criterion to be met. As mentioned in Ref. [19], the engine acceleration 
(fan/compressor speed changes) response to fuel flow characteristics can vary over a wide 
range due to the cumulative effect of various variables such as wear and differences in 
manufacturing tolerances, variable geometry setting, inlet distortion, environmental 
conditions etc. For the purposes of robustness analysis, this variation in propulsion 
dynamics can be represented as uncertainty in the A matrix and B matrix elements 
corresponding to the change in rotor speed rates (ft2, N25) with change in rotor speeds and 
fuel flow respectively, i.e. A(i,j) and B(i,k) with i and j = 9,10 and k = 7 in the design 
model of Eq. (1). Structured singular value robustness with the 80 Knot design model and 
the current controller design was performed considering independent multiplicative 
perturbations in these 6 elements of the A and B matrices. This analysis gave \i = 1.276 
which implies a stability margin SM = i = 0.79, i.e. the closed— loop system will remain 
stable for changes of up to 79 % from the 80 Knot model in these A and B matrix elements. 
Thus the control design provides robustness to large variations in the engine dynamic 
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response characteristics. The actual change in these elements between the design model and 
the off— design 60 Knot and 100 Knot models was much less than 79 % thus indicating that 
the design controller will provide robustness to changes in engine dynamics due to changes 
in the engine operating point within the transition envelope. 

Controller Scheduling 

The general approach to controller scheduling is to design linear controllers at 
various operating points and then perform some kind of curve fit to the various controller 
gains with the critical operating point conditions as the independent parameters. For a 
controller of order n with m inputs and r outputs, there will be n(l+m+r)+mr controller 
parameters to schedule with the controller represented in a minimal parameter state— space 
form [20]. For the current control design, n=14, m=18 (y of Eq. (4), 6 tracking errors and 
WB3) and r=12 resulting in 434 parameters to be scheduled. Even with further controller 
simplification such as removing the insignificant feedback paths, and partitioning the 
controller into separate airframe and propulsion sub— system controllers [1] with further 
partitioning of the airframe controller into longitudinal and lateral subcontrollers, the 
number of parameters to be scheduled may still be quite large. Moreover, the whole 
procedure of designing simplified controllers for the different operating points will be quite 
cumbersome and time consuming. A simpler controller scheduling scheme that exploits the 
robustness properties of the nominal controller is discussed in the following. 

The fundamental objective in using modern robust multivariable control design 
techniques is to reduce the controller complexity while guaranteeing the desired 
performance and stability robustness characteristics. A robust control design problem 
was formulated in this paper. The stability robustness of the resulting control design was 
demonstrated with special attention being given to variations in the elements of the system 
A matrix with changes in the operating point — trim airspeed for the model variations 
being considered. For a control design which is robust for changes in the plant A matrix, 
controller scheduling should only be needed to account for changes in the plant B matrix, 
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i.e. changes in control effectiveness. For example, for the aircraft models being considered, 
the effectiveness of the aerodynamic control surfaces decreases with airspeed. Therefore the 
control of the aircraft attitude angles and angular rates is transferred from aerodynamic 
surfaces to propulsion system generated moments via RCS and nozzle vectoring controls as 
the aircraft slows down. The desired control redistribution can be achieved by control 
scheduling of the form 

K(s) = K s K°(s) (8) 

where K(s) is the scheduled controller, K g is the control redistribution matrix and K°(s) is 
the controller designed for the nominal design point. Note that such a controller scheduling 
will remove the burden of having to do a nominal control design at too many operating 
points and also the number of parameters to be scheduled will be lower. For a controller 
with m outputs, the number of parameters to be scheduled will be at most m 2 which 
corresponds to 144 parameters for the current control design. However, the number of these 
parameters can be drastically reduced by taking into account the separation of the control 
effectors into different control axes. 

A procedure for developing controller scheduling of the form of Eq. (8) is discussed 
in the following and evaluation results are presented for the current design example. This 
controller scheduling approach was motivated by the control selector design approach of 
Ref. [21] wherein a decentralized, hierarchical approach to integrated flight propulsion 
control design was developed. 

The controller scheduling concept is demonstrated in Fig. 9 which shows the 
nominal (design point) control loop and the control loop at an off— design operating point. 
The plant feedthrough matrix D is not shown in Fig. 9 for simplicity of discussion. If the 
design controller is robust to changes in the A matrix and the controlled and feedback 
outputs, z and y respectively, are the same for off-design models, then choosing Kg such 
that 

BKg s B° 


( 9 ) 
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will ensure that the loop at point (i) in Fig. 9 is invariant to changes in B matrix thus 
guaranteeing robustness to changes in control effectiveness. One can then solve for an 
approximate value of Kg by using pseudo— inverse techniques [22] on Eq. (9). In order to 
perform the pseudo— inverse, it is important to ensure that the control effectiveness matrix 
B used in Eq. (9) is full rank, i.e. it represents the minimum number of independent 
controls for the problem being considered. For the integrated flight propulsion model of (1) 
to (3) the control degrees of freedom are the generation of the forces and moments for 
aircraft translational and rotational motion along the three axes (u,v,w,p,q,r), and the 
engine fan acceleration (N2). The controller scheduling gains are then given by 

Kg = N(EN)#E° (10) 

where E consists of the first six and the 9th row of B of Eq. (1) with all the columns, # 
represents pseudo— inverse, and N = diag(u ) with u corresponding to the maximum 

max max 

allowable control deflections. 

For the current design the controller scheduling gains Kg were obtained for the 
off— design 60 Knot and 100 Knot models using the appropriate E matrices. Fig. 10 shows a 
comparison of the closed— loop system response to lateral (Pv) command for the three flight 
conditions with the scheduled controller. Note that the controller scheduling maintains roll 
tracking performance with automatic turn coordination demonstrated by the low sideslip 
response. This controller scheduling provided robust decoupled command tracking for all 
the controlled variables in z of Eq. (5) except for the pitch rate variable Qv. For the case of 
Qv command, although the tracking performance with the scheduled controller was much 
improved over that without any controller scheduling, there was a large degradation in 
performance at the off— design conditions when compared to the nominal performance. To 
understand the reasons for this degradation in pitch tracking performance, further 
evaluations were performed for the Qv command case. First the closed— loop system was 
evaluated with the nominal controller keeping the plant B matrix constant at the nominal 
design condition and changing the A matrix to the off— design conditions. Then the 
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closed-loop system was evaluated with the scheduled controller keeping the A matrix 
constant at the nominal design condition and changing the B matrix to the off— design 
conditions. Robust tracking performance was obtained for the latter case but not for the 
former case. This result indicates that the controller scheduling scheme does provide 
robustness to changes in control effectiveness, however the nominal control design itself is 
not robust in pitch tracking performance to changes in the A matrix. This lack of 
performance robustness is due to the fact that aircraft pitch dynamics are a strong function 
of airspeed. Note that the results presented in the previous section had mainly 
demonstrated stability robustness and not performance robustness to changes in the A 
matrix. Clearly then guaranteeing performance robustness to changes in the A matrix is 
very critical to making the suggested controller scheduling scheme work. Although 
performance robustness analysis can be performed using structured singular value concepts, 
reliable numerical algorithms to do such analysis are not currently available. Since the 
internal noise model procedure was successful in building— in performance robustness in the 
lateral-directional axes, the possibility of using a similar approach for the pitch axis will be 
investigated in future design studies. 

Conclusions 

The Integrated Flight/Propulsion Control system design presented in this paper 

demonstrates the applicability of an control synthesis technique to integrated control 

design for complex systems such as Short Take-Off and Vertical Landing (STOVL) 

aircraft. The capability to address closed-loop performance, robustness and control 

actuation trade-offs within the framework of an H control problem formulation was 

00 

demonstrated. In particular, the control design for the IFPC problem, as formulated 
herein, provided robustness to the absolute nonlinearity associated with the compressor 
bleed flow demand from the Reaction Control System and to changes in the plant system 
dynamics over the transition flight envelope. A controller scheduling approach which has 
the potential to considerably reduce the complexity of control implementation over the 
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flight envelope and which exploits the robustness properties of the linear point design 
based controller was developed and demonstrated. 
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Table 1. Trim Control for Linear Models 


(Listed 

as percent (%) of 

80 Knot Model trim values) 

Vq (Knots) 

60 

80 

100 

A8 0 (%) 

11 

100 

268 

ETA 0 (%) 

96 

100 

118 

A78 0 (%) 

133 

100 

32 

WF q (%) 

116 

100 

71 

N2 0 (%) 

104 

100 

90 


2. Eigenvalues of Linear Models 



Eigenvalue 

60 

80 

— 7.9e-02 

— 8.5e — 02 

— 0.09±j0.24 

— 0.11±j0.28 

1.05 

1.29 

-1.64 

-1.73 

-1.74 

-2.09 

0.05±jl.99 

— 0.23±j2.27 

— 5.58±j0.74 

-4.12 


-7.11 

-29.73 

-29.39 

-38.67 

-38.21 

-172.0 

-199.7 



Descrintion 

100 

V Q (Knots) 

— 7.9e— 02 

Spiral Mode 

— 0.12±j0.31 

Phugoid Mode 

1.44 

Unstable SP a 

-1.77 

Roll Mode 

-2.45 

Stable SP 

— 0.48±j2.53 

Dutch Roll 

-3.33 

Rotor Speeds 

-3.83 


-27.73 

Temperatures 

-39.95 


-277.1 

Pressure Mode 


a Short Period 



Magnitude 


W c (s) 



Fig. 3 Block Diagram for Formulation of IFPC Feedback Design 
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Fig. 4 Singular Values of Scaled Design Plant 
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Fig. 5 Closed-Loop Weighted Norms for Controller 
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Fig. 6 Closed— Loop Tracking Singular Values with Full and Reduced— Order Controllers 
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Fig. 7 Pitch Trackin, 
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Fig. 8 Stability Robustness Measure p(M) for Variations in Airframe portion of Plant A 
Matrix 



(b) Off-Design Control Loop 


Fig. 9 Controller Scheduling for Variation in Control Effectiveness 
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Fig. 10 Roll Variable fPv) Tracking and Sideslip Response at the three operating points 
with Scheduled Controller 
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